Supercurrent and Andreev bound state dynamics in superconducting quantum point 

contacts under microwave irradiation 



F.S. Bergeret^ 2 ,P. Virtanen 3 - 4 , A. Ozaeta 1 , T.T. Heikkila 4 , J.C. Cuevas 51 

1 Centro de Fisica de Materiales (CFM-MPC), Centra Mixta CSIC-UPV/EHU, 
Manuel de Lardizbal 5, E-20018 San Sebastian, Spain 

2 Donostia International Physics Center (DIPC), 
Manuel de Lardizbal 4, E-20018 San Sebastian, Spain 

3 Institute for Theoretical Physics and Astrophysics, 
University of Wurzburg, D-97074 Wiirzburg, Germany 
i Low Temperature Laboratory, Aalto University, P.O. Box 15100, FI-00076 AALTO, Finland 
Departamento de Fisica Teorica de la Materia Condensada, 
Universidad Autonoma de Madrid, E-28049 Madrid, Spain 
(Dated: January 20, 2013) 

We present here an extensive theoretical analysis of the supercurrent of a superconducting point 
contact of arbitrary transparency in the presence of a microwave field. Our study is mainly based on 
two different approaches: a two-level model that describes the dynamics of the Andreev bound states 
in these systems and a fully microscopic method based on the Keldysh-Green function technique. 
This combination provides both a deep insight into the physics of irradiated Josephson junctions 
and quantitative predictions for arbitrary range of parameters. The main predictions of our analysis 
are: (i) for weak fields and low temperatures, the microwaves can induce transitions between the 
Andreev states leading to a large suppression of the supercurrent at certain values of the phase, 
(ii) at strong fields, the current-phase relation is strongly distorted and the corresponding critical 
current does not follow a simple Bessel-function-like behavior, and (iii) at finite temperature, the 
microwave field can enhance the critical current by means of transitions connecting the continuum 
of states outside the gap region and the Andreev states inside the gap. Our study is of relevance 
for a large variety of superconducting weak links as well as for the proposals of using the Andreev 
bound states of a point contact for quantum computing applications. 

PACS numbers: 74.40.Gh, 74.50.+r, 74.78.Na 



I. INTRODUCTION 



In 1962 Josephson predicted that a dissipationless cur- 
rent (supercurrent) could flow in a junction between 
two superconductors (S) weakly coupled by an insulat- 
ing barrier fi which was confirmed experimentally shortly 
afterwards by Anderson and Rowell.- Soon after this con- 
firmation, it became clear that this phenomenon, referred 
to as the dc Josephson effect, could take place in a vari- 
ety of superconducting weak links such as Dayem bridges, 
SNS junctions, where N corresponds to a normal metal 
bridge, or large point contacts^^ The only difference be- 
tween these systems lies on the exact current phase rela- 
tion (CPR), which depends on the characteristics of the 
constriction linking the superconducting leads.— 

In recent years, the dc Josephson effect has 
been investigated in novel superconducting junc- 
tions with weak links based on atomic contacts^— 
carbon nanotubes^— fullerenes^ semiconductor 
nanowireS ) 13 i 14 or grapheneJ^— Some of these nanos- 
tructures fall into the category of a superconducting 
quantum point contact (SQPC), where the constriction 
has a length much smaller than the superconducting 
coherence length. In this limit, and in the absence of 
strong interactions in the constriction, the dc Josephson 
effect can be described in a unified manner using 
two basic concepts of mesoscopic physics, namely the 



concepts of conduction channels and Andreev bound 
states. In the normal state, the coherent transport 
through a mesoscopic system can be described in terms 
of the independent contributions of the eigenfunctions 
of the transmission matrix from the structure, known 
as conduction channels, and these contributions are 
determined by the corresponding transmission coeffi- 
cients {ti}. In the superconducting state, the electrons 
(holes) transmitted in a conduction channel are Andreev 
reflected at the electrodes as holes (electrons) in the 
same channel. This process is successively repeated in 
both electrodes leading to the formation of a pair of 
bound states in the gap region. These are known as 
the Andreev bound states (ABSs). In the case of a 
single-channel SQPC with transmission r, the energies 
of the ABSs are given b y 18 ' 19 



Ef(<p,T) = ±E A (<p,T) = ±Ay/l - rsin%/2), (1) 

where A is the superconducting gap and ip is the phase 
difference between the order parameters on both sides 
of the junction. In equilibrium, these two states carry 
opposite supercurrents I^(<p) — {2e/h)dE^/dip, which 
are weighted by the occupation of the ABSs (determined 
by the Fermi function). In the case of a multichannel 
SQPC, the supercurrent is simply given by the sum of 
the contributions from the individual channels^ 9 - 

This unified microscopic picture of the dc Josephson 
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effect has been confirmed experimentally in the con- 
text of atomic contacts by Delia Rocca and coworkers, 8 
In particular, these authors measured the CPR of an 
atomic contact placed along with a tunnel junction in a 
small superconducting loop and found an excellent agree- 
ment with the theory using the independently determined 
transmission coefficients. At this stage, one may wonder 
whether it is possible to control the occupation of ABSs 
of a SQPC with an external field, and in turn to control 
the supercurrent. This is the main issue explored in this 
work and for this purpose, we present here an extensive 
theoretical analysis of the supercurrent and the dynam- 
ics of the ABSs of a SQPC under microwave irradiation. 
This is a basic problem in mesoscopic superconductivity, 
which is also relevant for the field of quantum comput- 
ing since the ABSs of a SQPC have been proposed to be 
used as the two states of a qubit.— ~— In this proposal, a 
microwave field can be used for the spectroscopy of the 
two-level system or to probe its quantum state by current 
measurements. 

The microwave- assisted supercurrent in SQPCs is of- 
ten discussed in the framework of the adiabatic approxi- 
mation (see Section II), where one assumes that the ABSs 
follow adiabatically the microwave field. This approxima- 
tion does not take into account the possible transitions 
between the ABSs and therefore, if fails to describe the 
current at high frequencies or for highly transmissive con- 
tacts, where the energy difference between the states can 
be rather small. The first microscopic analysis of this 
problem for a SQPC of an arbitrary transparency was 
reported by Shumeiko and coworkers^, These authors 
studied the limit of weak fields and predicted the pos- 
sibility to have a large suppression of the current due 
to resonant transitions between the ABSs. Later, other 
aspects of this problem, including the dynamics of the 
ABSs, have been addressed focusing on the linear re- 
sponse regime^ - — A complete solution of this problem, 
valid for an arbitrary range of parameters, has only been 
reported very recently^ In this latter work, we devel- 
oped a theory of the supercurrent through a microwave- 
irradiated SQPC in the framework of the Keldysh-Green 
function technique. This theory allowed us to put for- 
ward new predictions such as the evolution of the CPR 
with the radiation power and the possibility to enhance 
the critical current at finite temperatures by irradiating 
the junction. Here, we describe in detail this theory (see 
Section IV) and, in particular, we present new analyt- 
ical results that elucidate the origin of the microwave- 
enhanced supercurrents. 

In the process of understanding the results of the exact 
theory, we are confronted with the question of to what 
extent the physics of microwave- irradiated SQPCs can be 
understood in terms of just the dynamics of the ABSs, 
i.e., in terms of a natural extension of the argument de- 
scribed in the previous paragraphs for the case of a junc- 
tion in equilibrium. In order to answer this question, we 
make use of the two-level Hamiltonians of a SQPC exist- 
ing in the literatur o 21 ' 28 and we compare the results with 



the exact theory. This comparison serves in turn to es- 
tablish the range of validity of these two-level models. It 
is worth stressing that within these models the computa- 
tion of the dc properties such as the supercurrent or the 
average occupation of the ABSs for arbitrary radiation 
is a highly non-trivial task. In order to carry it out, we 
have developed a new powerful method which allows us 
to compute any dc quantity in an arbitrary two-level sys- 
tem driven out of equilibrium by a periodic perturbation. 
This method is described in Section III and it constitutes 
one of the main results of this work. With the help of this 
method, we show that with the Hamiltonian of Ref. Hi] 
one can nicely reproduce the exact results at low temper- 
atures and low radiation powers. Moreover, this analysis 
allows us to obtain analytical results for the supercurrent 
dips produced by microwave-induced transitions between 
the ABSs. 

The rest of the paper is organized as follows. In the 
next section we briefly review the equilibrium properties 
of a SQPC as well as the basic results of the adiabatic 
approximation. In Section III we study the dynamics 
of the ABSs under a microwave field within the two-level 
Hamiltonian of Ref. [U In particular, we describe a novel 
method that allows us to obtain the CPR for any power 
and frequency of the external field, and we also derive 
analytical expressions for the supercurrent beyond the 
rotating wave approximation. In Section IV we discuss 
the Keldysh-Green function technique, which describes 
the supercurrent for an arbitrary range of parameters, 
including also the contribution of the continuum of states 
outside the gap region. We present a detailed comparison 
of the results of this technique at zero temperature with 
those obtained with the two-level model. Moreover, we 
analyze in detail the phenomenon of microwave-enhanced 
supercurrent at finite temperatures, for which we present 
analytical results. Finally, Section V is devoted to some 
additional discussions and to summarizing the main re- 
sults of this work. 



II. SYSTEM AND ADIABATIC 
APPROXIMATION 

We consider a SQPC consisting of two identical super- 
conducting electrodes with an energy gap A, linked by a 
single conduction channel of transmission r. Our main 
goal is to compute the supercurrent through this system 
when it is subjected to a monochromatic microwave field 
of frequency u). We assume that the external radiation 
generates a time-dependent voltage V{t) — Vosmivtr- 
where the amplitude Vq depends on the power of the 
external radiation source, and eventually also on the po- 
larization of the radiation. According to the Josephson 
relation, this voltage induces a time-dependent supercon- 
ducting phase difference given by 

4>(t) = (p + 2a cos tut, (2) 
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where ip is the dc part of the phase and a = eVb /hui is 
parameter that measures the strength of the coupling 
the electromagnetic field and it is used here as a param 
ter to be determined by comparing with the experiment 
As explained in the introduction, in the absence of n 
crowaves the supercurrent can be expressed as a su 
of the contributions of the two ABSs as I e q(<p) 
I^Hf(E^) + I^np(Ej L ), np(E) being the Fermi disti 
bution function, which yield: 



,20 



Ieq{<P) 



eA 2 rsiny ^ / E A (ip) \ 
2h E A (<p) tan±1 1, 2k B T J 



( 



where Ex is defined in Eq. (fTJ) and T is the temper 
ture. In the tunnel regime (r <C 1), this expression r 
duces to the sinusoidal CPR given by the Ambegaoka 
Baratoff formula: 30 I C q{f) = Ic sm <Pi with Iq 
(eAr/2fi) tanh(A/2fcBT). At perfect transparency (r 
1), this expression reproduces the Kulik-Omelyanchi 
formula*^ I C q(<fi) = h sin(<p/2) tanh(Acos((^/2)/2fc B T 
Here, Iq = eA/h is the zero-temperature critical curre: 
for t = 1 and we frequently use it below to normali 
the supercurrent in the different graphs. According 
Eq. (|3|), at zero temperature only the lower ABS co 
tributes to the supercurrent I oq = it, while at a fini 
temperature the negative contribution from the upp 
ABS leads to a decrease of the total supercurrent. 

The simplest approach to compute the supercurrent 
the presence of the microwave field is the so-called adi 
batic approximation. In this approximation one assumes 
that the ABSs follow adiabatically the ac drive and there 
are no direct transitions between them. Thus, the CPR 
in this approximation is obtained by replacing the sta- 
tionary phase ip in Eq. (J3|) by the time-dependent phase 
4>(t) of Eq. (|2|) , which leads to the following result 



oo 



Iad((p, a) = ^ InJo(2na) sin(n^), 

n=l 



(4) 
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FIG. 1. (Color online) Panels (a) and (b): The current-phase 
relation in the adiabatic approximation for r = 0.2 (a) and 
t — 0.95 (b). The different curves correspond to different 
values of a as indicated in the graphs. The current is given in 
units of To = eA/h, where A is the value of the superconduct- 
ing gap at T = 0. (c) The zero-temperature critical current as 
a function of a for three different values of the transmission 
r. Notice that critical current is normalized by its value in 
the absence of microwaves. 



where /„ = (1/V) f^* dtp I cq (<p) sm(ntp) are the harmon- 
ics of the equilibrium CPR of Eq. ([3]) and Jo is the zero- 
order Bessel function of the first kind. Notice that the 
current in this approximation does not depend explic- 
itly on the radiation frequency. We illustrate the results 
of this approximation in Fig. [1] for the zero-temperature 
case. In particular, in the two upper panels we show 
the CPR (obtained from Eq. Q) for two different trans- 
missions and several values of the a parameter (related 
to the microwave power). Panel (a) corresponds to the 
tunnel limit (r = 0.2) where the CPR is sinusoidal ir- 
respective of the radiation power, while in panel (b) we 
show the results for a high transmission of t = 0.95. In 
this latter case, the critical current is reached at different 
values of the phase depending on the value of a. Notice 
that no matter the value of the phase <p, the magnitude of 
the supercurrent is always suppressed by the microwaves 
as compared with the zero-field result (a = 0), which is 
true at any temperature. With respect to the behavior 



of the critical current Ic(a), as one can see in Fig. [He), 
it decays in a non-monotonic manner, as governed by the 
Bessel function Jq. 



III. THE TWO-LEVEL MODEL 

It is instructive to start our analysis towards a mi- 
croscopic theory by restricting ourselves to the study 
of the contribution of ABSs, ignoring for the moment 
the continuum part of the spectrum. This can be done 
with the help of the two-level models that have been de- 
rived in Refs. and [2l| to describe the dynamics of a 
SQPC under external ac fields. The models of these 
two references coincide at equilibrium, but they differ 
slightly when the phase depends on time. In particular, 
the model of Ref. [2l| ensures charge neutrality, while the 
model of Ref. [28| does not. For this reason, we base our 
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discussion here on the model put forward by Zazunov 
and coworkers^ In this model, the SQPC is described 
by the 2x2 Hamiltonian 



H B (t) = Ae~ i& ^ 2 f cos Za, + r sin |<r y 



(5) 



where r = y/1 — r and <p(t) is the time-dependent phase 
given by Eq. (J2J. This Hamiltonian is written in the 
ballistic basis of right- and left moving electrons, which 
are eigenvectors of the current operator in the perfectly 
transmitting case (r = 1). For our subsequent analy- 
sis it is more convenient to work in the instantaneous 
Andreev basis {|+)0(t)j |— )<p(t)}, whose basis vectors are 
time-dependent. This is the basis where the Hamiltonian 
of Eq. ([5]) becomes diagonal in equilibrium. The Andreev 
basis is obtained from the ballistic basis by means of a 
transformation H A {t) = R) (t)H B {t)R(t) generated by 
the unitary matrix 



R(t) = e -i*xr$ e -i%&* e -i9{<l>)& y 



(6) 



where 6(<fi) = (l/2)arctan[rtan(</>/2)]. With this trans- 
formation the Schrodinger equation for a state vector 
= (a{t),f3(t)) T becomes 



where 

H A (t) = E A (ct>(t))a z - 



i&tf (t) = H A (*)*(*) , 

rrA 2 sin 2 (0(t)/2) 



4[£a(0W)] 



(7) 



<p(t)a v , (8) 



and (}>(t) — d<p(i)/dt. Moreover, in the previous two 
equations, and in the rest of this section, we set % = 1. 
The corresponding current operator can be written as 



I A (t) = 2eE' A {(j>{t))a z 



err A 2 sin 2 (0(t)/2) „ 



(9) 



where the prime in E' A means derivative with respect to 
the argument (the time-dependent phase in this case) . To 
obtain the expectation value of the current at different 
times, Eq. ([7]) needs to be solved. Despite the apparent 
simplicity, this task has nontrivial aspects: straightfor- 
ward numerical approaches run into problems, as both 
very fast (t^ 1 ~ w) and very slow (t^ 1 ~ E A — nuj) 
time scales can be simultaneously present. No closed- 
form analytical solution can be obtained either^ 2 - and 
the significantly nonlinear coupling to the drive makes 
it more difficult to derive approximations via standard 
routes 

Focusing our analysis on time-averaged quantities, we 
can obtain accurate analytical and numerical results via 
a systematic Floquet-type approach. We are interested 
in two physical quantities: the dc current 



I = lim — 

t— too t 



(10) 



and the time-averaged populations of the Andreev levels 



P± = lim 7 



1±<7 2 



*(*0 • (ii) 



Below, we show how to obtain 7, although the method 
described can as well be used to compute any other time- 
averaged quantity, including p± . 

We first introduce a modified Hamiltonian 



H A (t, X ) = H A (t)+ X iA(t) 



(12) 



where x is a parameter conjugate to the observable, and 
it is set to zero at the end of the calculation. The solution 
of the Schrodinger equation ^(t, \) can be formally writ- 
ten by introducing the time evolution operator J7(i, 0; x) 

*(t, x) = Te- 1 /o dt'H A {t', x )y Q = jy^ . x) ^ q ( (13) 

where T indicates time ordering and is the state vec- 
tor at t — 0. We define now the generating function: 

S(t, X ) = *tU(0, t; X = 0)U(t, 0; X )*o . (14) 

One can easily check that the dc current defined in 
Eq. (fTOj) can be written as 



I 



lim ~d x S(t,x)\ x =o 

t— >oo t 



(15) 



Thus, we only need to compute the function S, or, equiv- 
alent^, the evolution operator U(t, 0;x) = U(t;x)- 

Since our Hamiltonian is periodic in time with a period 
T = 2it/uj, i.e., H A (t,x) — H A (t + T;x), we can define 
two periodic (Floquet) states v± via the eigenvalue prob- 
lem 



^(T;x)v ± (x) = e ±lB ^ T v ± (x). 



(16) 



The symmetry of the two eigenvalues follows here from 
the fact that H A (t, x) and log[[/(T; x)] are traceless, and 
U(T;x) is unitary. Moreover, from the periodicity of the 
Hamiltonian it follows that 



U(nT; X ) = U(T;xT = V( X )e 



iE( X )nTa^y-l 



V- L (x), (17) 



where the eigenvectors v± form the columns of the uni- 
tary matrix V. Replacing in Eqs. (|T4^) and (TTSl) t by nT 
and taking the limit n — > oo we now find the derivative 
with respect to x : 



1 

nT 



d x U(nT-x) 



iV- x {x)OzZ lE{x)nT&z V(x) 



dE{ X ) 
dx 



Thus, the dc current is given by 



I = 



-*S( 



dE( X ) 



(18) 



(19) 



This exact expression for the dc current is very useful for 
numerics. It is easy to compute and it handles the fast 
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and slow time scales of the problem separately. In order 
to obtain the dc current, one needs first to integrate the 
Schrodinger equation with the Hamiltonian of Eq. ([12"]) 
over one period to find the 2x2 matrix U(T; x) , then one 
computes its eigenvalues ±E and eigenvectors v±, and 
finally the derivative d x E(x) is computed via numerical 
differentiation. 

In order to have a first impression of the results from 
this two-level model, we show in Fig. [2] a few examples 
of the CPR of a highly transmissive channel (r = 0.95) 
computed with the numerical recipe that we have just 
described^ The upper panels of this figure show the 
CPR for a moderate power (a = 0.15) and three different 
values of the microwave frequency. For comparison, we 
also show the result obtained with the adiabatic approx- 
imation of Eq. (0}. As one can see, the main difference 
is the appearance in the results of the two-level model of 
a series of dips at certain values of the phases where the 
current is largely suppressed. It is easy to understand 
that such dips are due to microwave-induced transitions 
between the ABSs. These transitions enhance the pop- 
ulation of the upper ABS, which at zero temperature 
would be empty otherwise, and at the same time they 
reduce the occupation of the lower ABS. This redistribu- 
tion of the quasiparticles in turn results in a suppression 
of the current. The microwave-induced transitions occur 
with the highest probability when the distance in energy 
between the ABSs (the Andreev gap) is equal to a multi- 
ple of the photon energy, i.e., when 2Ea{<p) = nuj, where 
n = 1,2,... is the number of photons involved in the 
transition. If this condition is expressed in terms of the 
phase ip, it adopts the from 



tonian 



obeys 



<Pr. 



2arcsinv / [l - (nw/2A) 2 ]/r, n=l,2,.. 



(20) 



A detailed analysis shows that this expression reproduces 
the positions of all the dips appearing in the examples of 

Fig.m 

This interpretation of the origin of the dips in the CPR 
can be corroborated by a direct analysis of the occupa- 
tions of the ABSs. Following the same numerical recipe, 
we have also computed the average occupation of the up- 
per ABS, p+j for the examples shown in the upper panels 
of Fig. O The results can be seen in the lower panels of 
this figure and, as one can observe, there is a clear one- 
to-one correspondence between the current dips and the 
enhancement of the time-averaged population of the up- 
per state. In particular, whenever the upper state reaches 
a population equal to 1/2, the current vanishes exactly. 

The method described above is not only very conve- 
nient for numerical calculations, but it also provides a 
route to obtain analytical results. In what follows, we 
show how this method can be used, in particular, to gain 
a further insight into the microwave-induced supercur- 
rent dips. To proceed, it is useful to first rewrite Eq. ([19]) 
in a more convenient form. In particular, we would like 
to avoid the calculation of eigenvectors in this equation. 
This can be done by noting that the unperturbed Hamil- 



a x Hp^& x = -Ha ■ 



(21) 



Consequently, v_ cx a x Y + , and the dc current given by 
Eq. ([19"]) can be written as 



/ = -vi (p - a x p a x ) v + 



dE( X ) 



dx 



(22) 



where po = ^o^o- Using the expression for the change of 
an eigenvalue due to a perturbation, we can finally write 



dE( x ,p) dE(x,n) 



dp, 



dx 



Here E(x> A*) is an eigenvalue of the matrix 
M (x, M) = j, M{T, x) + V (A) - VxPoVx) 



(23) 



(24) 



0(T, x) = log[Z7(T, x)] and p is an additional perturba- 
tion parameter. The problem is now reduced to finding 
the eigenvalues of a 2x2 matrix. 

As discussed above, the dc current for weak fields only 
deviates from the adiabatic result close to the resonant 
conditions nu — 2Ea (with n = 1,2,3,...), where the 
transitions between the ABSs are more likely. In order 
to study what happens close to these resonant situations, 
we can consider the problem in a rotating frame, and 
rewrite the evolution operator U(t) defined in Eq. ([13]) 
as 

U(t) = e -iW n t Te -if*dt'H n (t') = e -^n*u n (t) , (25) 



where W n 
is 



moa z /2 and the rotating- frame Hamiltonian 



H n {t) = 



„iW n t 



H A (t) - W n ]e 



-iW n t 



(26) 



The generating function can then be written as in 

Eq. (fl"4"]) simply by substituting H\ by H n . The ad- 
ditional exponential factors simply cancel out, and the 

Hamiltonian H n (t) remains periodic. Thus, we can pro- 
ceed exactly as above. 

The key idea that allows us to obtain analytical results 
is the fact that for weak fields (a -C 1), the dynamics in 

the rotating frame are slow (H n is small) around the cor- 
responding resonance. For this reason, we can use the 
Magnus expansion 3 — to determine the matrix Cl appear- 
ing in Eq. 



n(T) 



dhH n (h) 



(27) 



dU 



H n (ti), H n (ti) 



This is essentially an expansion in the parameter A„ ~ 
2wk(Ea — nuj/2)/oj, which indeed is small close to a res- 
onance. 
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FIG. 2. (Color online) (a-c) Zero temperature supercurrent, in units of Io = eA/h, as a function of the phase for r = 0.95, 
a — 0.15 and three different values of the microwave frequency, as indicated in the upper part of the graphs. The solid lines 
correspond to the numerical results obtained with the two-level model, while the dashed lines is the result obtained with the 
adiabatic approximation, (d-f ) Time-average occupation of the upper ABS for the cases shown in the upper panels. 



We proceed now computing the dc current close to 
the first resonance u> = 2E\, assuming that initially the 
system is in its ground state VfJ = (0, 1). We choose W = 
w(T z /2, and take only the first term of the expansion (|27J). 
expanding up to the first order in a and x- The time 
integral is straightforward to evaluate, and we obtain 



lil 1 (T,x)^[E A -~+2exE' A \a z 



(28) 



-a[(A 2 -E 2 A )^-x(A 2 + E 2 A )2eE' A ]a x 



Note that this expression is analogous with the well- 
known rotating wave approximation, with the differ- 
ence that by considering the generating function, our 
formalism takes the time dependence of the operator 
I(t) into account. For the eigenvalues of the matrix 

Mi = (i/T)Qi + fia z we obtain 



E 2 = [E A ~-+n + 2e X E' A ] 2 



(29) 



■r (ra/2^) a [(A a - E\)^ - X (A 2 + E 2 A )2eE' A ] 2 . 

Finally, working in the limit (w — 2E A )/A <C 1 for sim- 
plicity, we find the dc current from Eq. 



h{y,u,cc)tt-2eE' A 1 



r.l 



(30) 



where the resonant frequency uj\ = 2E A equals the un- 
perturbed Andreev level spacing 2E A (up to first order in 



a), and f2 r ,i = raoj(A 2 — E\)/2E\ is the corresponding 
Rabi frequency. This expression tells us that the current 
vanishes exactly at the resonant condition u) = u)\ and 
that the width of the current dip is given by f2 r .i, which 
is linear in a. Moreover, its form clearly suggests that the 
populations of the two states undergo Rabi oscillations 
with the frequency f2 r ,i, as usual in two-state systems, 
and the time-averaged populations of the ABS coincide 
at the resonance. As a consequence, the dc current drops 
to zero at the resonance, a result that qualitatively coin- 
cides with the prediction in Ref. [2^. 

We can also determine the dc current at the higher res- 
onances, for example for w Ea- In this case we work 
in the frame corresponding to W2 — w<r z - As the reso- 
nance is due to two-photon processes, terms up to order 
a 2 must be taken into account, which requires including 
the first two terms in Eq. (|2"T|) . The computations are 
again straightforward, and up to the second order in a 
we obtain 



lfi a w [E A - + a 2 E" A + 2e X E' A + r 2 a 2 u ^j^- ] a z 



(31) 



ra 2 E' A u 



2E\ 



•E A 



l)-E A )]a x 



For simplicity, we dropped terms of order ax, which do 
not essentially affect the form of the resonance. As above, 

the current is obtained from the eigenvalues of M2 = 
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done in Fig. [3] where we have considered the case of a 
weak microwave field (a = 0.05). There is an excellent 
agreement between Eq. (|33|) and the numerical results, 
apart from the fact that the numerics also include a dip 
produced by three-photon processes, which we have left 
out from the above approximation. 

We can conclude this section by saying that in spite 
of the simplicity of the two-level model considered here, 
such a model captures the essential physics of the 
microwave-irradiated SQPC and, as we show in the next 
section, it provides accurate results for not too high fre- 
quencies and up to moderate radiation power. As we 
establish in the next section, the limitations of the two- 
level model are mainly related to the fact that it does 
not take into account the contribution of the continuum 
of states outside the gap region. 



IV. THE KELDYSH-GREEN FUNCTION 
APPROACH 



FIG. 3. (Color online) (a) Current-phase relation in the two- 
level model, obtained from Eq. (|33p (solid line) and numerics 
based on Eq. (|19|) (dotted line). Parameters are uj = 0.6A, 
r = 0.95, and a = 0.05. The n = 3 resonance is not included 
in the analytical approximation. Inset: close-up of the second 
resonance, (b) Time-averaged population p+ of the upper 
Andreev state for the same parameters. 



(i/T)Q,2 + H&z, and it adopts the form 
I 2 ((p,u},a) 



2eE A ( 1 - 7 "ff n2 

1 (w - u 2 y + n 2 



(32) 



where u 2 = £-a 
0,2 = ra 2 (2A 2 



- ol 2 E" a + r 2 a 2 (A 2 - E 2 A ) 2 /12E 3 A and 
E A )E' A /2E A . Here, one can observe 
that the resonant frequency is shifted from the position 
2uj = 2Ea by two contributions: the first arises from 
nonlinearities, and the second is the Bloch-Siegert shift. 33 
One can also go further and compute the dc cur- 
rent around resonances n > 2, although this gets pro- 
gressively more cumbersome as an increasing number of 
terms are required in the Magnus expansion, reflecting 
the increasing number of allowed multiphoton processes 
generated by the nonlinearities. One can however see 
from Eqs. (|3QI32p . and also check for higher resonances, 
that the width of the resonances scales with a™. More- 
over, one can show that within this model, the time- 
averaged current vanishes exactly at each resonance, i.e., 
I\ u = Un = 0. 

The results Eqs. (|3QI32[) can be combined into a single 
approximate expression 



-2eE' A 1 



O 2 



(w-wi) 2 + f7 2 1 



1- 



^,2 



(lu - lo 2 ) 2 + n 2 . 



( 33 ) 

The quality of this approximation can be established by 
comparing it with the exact numerical results. This is 



In the previous section we have analyzed the supercur- 
rent assuming that the only contribution comes from the 
ABSs. While this is true for a SQPC in equilibrium, it is 
not obvious that this should be the case in the presence 
of a microwave field. Indeed, at high frequencies or at 
high radiation powers, and especially at finite tempera- 
tures, transitions between the ABSs and the continuum 
of states outside the gap become possible and, in princi- 
ple, they can also contribute to the current. Therefore, in 
order to describe the complete phenomenology of irradi- 
ated SQPCs we must develop a fully microscopic theory. 
This is the goal of this section. 

Our microscopic theory is based on the Keldysh-Green 
function approach. In this approach the starting point 
is the expression for the quasiclassical Green functions 
of the left (L) and right (R) electrodes. In our case, 
these Green functions can be expressed in terms of the 
equilibrium Green functions g(t — t') as 



G\R(£)(M0 



_ p ±i4>(t)t 3 /2 



g{t - t')e 



=Fi0(t')r 3 /2 



(34) 



Here <j>(t) is the time-dependent phase given by Eq. ([2]) 
and the upper (lower) sign in the exponents corresponds 
to the R (L) electrode. The symbol " indicates that the 
Green functions are 4x4 matrices in the Keldysh-Nambu 
space, where they have the structure 



G = 



K 



G R G 
G 1 



(35) 



Here the symbol " indicates that the different elements are 
2x2 Nambu matrices. The retarded (R), advanced (A) 
and Keldysh (K) components of the equilibrium Green 
functions appearing in Eq. (|34[) are given by 



m = 



d ^e-^g{E), 

Z7T 



(36) 
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where 



gR(A) {E) = g R(A) {E)h + f R(A) {E)ih (37) 

g K (E) = [g R (E) - g A (E)} tanh(£/2fc B T) (38) 



and 
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E 



E 



y/(E ± ir\) 2 - A 2 A 



f R ( A \E), (39) 



where rj describes the inelastic scattering energy rate 
within the relaxation time approximation and T is the 
temperature. 

Different authors have shown that the transport 
properties of a point contact with an arbitrary time- 
dependent voltage can be described by making use of 
adequate boundary conditions for the full quasiclassical 
propagators These boundary conditions can be ex- 
pressed in terms of a current matrix 



jR j 

P 



K 



(40) 



which for the case of a single-channel SQPC of trans- 
mission r can be expressed in terms of the lead Green 
functions of Eq. (|34j) as2& 



I(t,t')=2r[G L ,G R ] o o[4-r(2-{G L ,G R } o )] 1 (t,t'). 

. (41) 

Here, the symbol o denotes the convolution over inter- 
mediate time arguments. Finally, the electric current is 
obtained by taking the trace 



I(t) = -Trr 3 I K (t,t) , 



(42) 



where T3 is the third Pauli matrix in Nambu space. 

Due to the periodic time dependence of the phase 
[Eq. @], the Green functions G L rm, and any products 
of them, admit the following Fourier expansion 

G(t,t')= e imut ' I ^-e-^-^^GomiE), (43) 



m— — 00 



where G nm (E) = G{E + nftw,E + mtuS) are the corre- 
sponding Fourier components in energy space, and n, m 
are integers. In particular, the Fourier components of 
G L rm can be deduced from Eq. (f3~4"|). For instance, for 
the left electrode, G nm (E) is given by 



(Gi)„ m — ^ r n igit* m , 
1 



where 
f 



r r 

f„ 



,r r , 



v. 



nm u 

VI 



(44) 



(45) 



Here, V nm = (i) m - n J m _„(a/2)e^/ 4 , where J n is the 
Bessel function of order n, and g n — g(E + nhu) is the 



equilibrium Green function matrix with the argument 
shifted in energy. 

From this discussion, it is easy to understand that the 
current adopts the general expression 



I(t) 



E 



■Ml 



(46) 



which means that the current oscillates in time with the 
microwave frequency and all its harmonics. These cur- 
rent components can be computed from the Fourier com- 
ponents in energy space of / in Eq. (|4T|). From that 
equation, it is straightforward to show that the Fourier 
components of I K are given by 



fK _\^,ARj>K,AKj>A 



(47) 



Here, we have defined the matrices A nm = 2t[Gl, GR\ nm 
and X nm = [41 - t(2 - {G L , Gr})}~^, which can be 
determined from the Fourier components of G_l(_r). Once 
the components of I K are obtained from Eq. (|47|) . one 
can compute the current. We are only interested here in 
the dc component, which reads 



I(tp,u,a) 



e 

4fi 



dE 



TrT 3 iQ (E,ip,uj,a). 



(48) 



The dc current can be calculated analytically in certain 
limiting cases: for example in the absence of microwaves, 
where it reduces to Eq. in the tunnel regime or for 
very weak fields. However, for arbitrary radiation power 
one needs to evaluate Eq. (j48|) numerically. In the next 
subsections we present the results for the dc current of 
this microscopic theory and we compare them with those 
obtained from the two- level model of Section HID 



Zero-temperature limit: Comparison with the 
two-level model 



We focus first on the analysis of the results of the exact 
theory at zero temperature. This allows us, in particular, 
to make a comparison with the two-level model of Section 
III and to establish its range of validity. 

In Fig. @] we show several examples of the CPR cal- 
culated with the microscopic approach (solid lines) for a 
highly conductive channel (r = 0.95) for several frequen- 
cies and low values of the radiation power (a <C 1). For 
comparison, we also show the results of both the two- level 
model (dashed lines) and the adiabatic approximation 
(dotted lines). As one can see, the main deviation from 
the adiabatic results is the appearance of a series of dips, 
as discussed in section III. These features, which origi- 
nate from the microwave-induced transitions between the 
ABSs, are accurately reproduced by the two-level model 
(both the position and the width of the dips) . There is a 
small discrepancy between the exact result and those of 
the two-level model for phases close to it, i.e., when the 
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FIG. 4. (Color online) Four examples of the zero-temperature 
current-phase relation for r = 0.95 obtained from the micro- 
scopic model (solid lines), the two level model (dashed lines) 
and the adiabatic approximation (dotted lines). The param- 
eters characterizing the microwave field are indicated in the 
different panels. 



level spacing between the ABS is very small. This is un- 
derstandable since the model assumes that tkj)(t) <§; 2Ea, 
which is not fulfilled when ip ~ it and r is close to 1. No- 
tice also that for the high-order dips (due to high-order 
photonic processes) , the suppression of the current in the 
two-level model is larger than in the case of the exact the- 
ory. The reason is the additional broadening introduced 
by the finite inelastic scattering rate used in the calcula- 
tions with the microscopic theory, which in this case is 

7] = icr 4 A. 

The good agreement between the microscopic theory 
and the two-level model in these examples can be under- 
stood as follows. At zero temperature, the lower ABS 
is fully occupied, while the upper one is empty. There- 
fore, for small values of a and hu < A transfer of quasi- 
particles between the continuum and the ABSs is not 
possible. The agreement between these models is fur- 
ther confirmed in Fig.[5ja), where the CPR is shown for 
hw = 0.6A, a = 0.1 and two lower values of the trans- 
mission (t = 0.6 and 0.8). In this case, the agreement is 
almost perfect for all phases. The reason is that now the 
smallest energy gap between the ABSs, which occurs at 
ip = n, is large enough to avoid the overlap of the levels in 
the presence of the microwave field. If the transmission 
is further reduced, no transitions can occur between the 
Andreev states and the adiabatic approximation becomes 
exact. 

From the discussion above, we can conclude that the 
two-level model provides an excellent description of the 
supercurrent at zero temperature and for weak fields 
(a -C 1). However, as the radiation power increases, 
the situation changes. This is illustrated in Fig. [SJb), 
where we show the CPR for a highly conductive channel 
(t = 0.95), a frequency haj = 0.3A and two values of 
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FIG. 5. (Color online) (a) The current-phase relation for a = 
0.1, ftw = 0.6A and two values of the transmission coefficient, 
r = 0.8 and r = 0.6. (b) The current-phase relation for 
huj — 0.3A, r = 0.95 and two values of a, 0.2 and 0.6. In both 
panels the solid lines correspond to the microscopic theory and 
the dashed lines to the two-level model. 



a. As one can see, the deviations between the results of 
the two-level model and the microscopic theory become 
more apparent as the power increases. The main reason 
for this discrepancy is the occurrence of multiphotonic 
processes, which become more probable as the power in- 
creases. These processes induce quasiparticle transitions 
between the ABSs and the continuum part of the energy 
spectrum, which are not included in the two-level model. 

As one could already see in Fig. EJb), as the radiation 
power increases the supercurrent dips broaden and the 
CPR acquires a very rich structure. We illustrate this 
fact in more detail in Fig. |5]where we show the evolution 
of the CPR with a for two values of the transmission 
(0.95 and 0.8) and for frequency fku = 0.3A. Notice that 
as the power increases, the dips disappear, the CPRs be- 
come highly non-sinusoidal, and in some regions of the 
phase the current can reverse its sign. These results are 
clearly at variance with those found within the adiabatic 
approximation (see Section II). They are a consequence of 
a complex interplay between the dynamics of the ABSs, 
which are broadened by the coupling to the microwaves, 
and the multiple transitions induced between the ABSs 
and the continuum of states. This very rich behavior 
has also important implications for the critical current, 
which for high transmission strongly deviates from the 
standard behavior described by the adiabatic approxi- 
mation. This is discussed below in detail. Finally, it is 
worth stressing that the values of a used in Fig. [5] are 
easily achievable in experiment, as demonstrated in the 
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FIG. 6. (Color online) The zero-temperature current-phase 
relation for hw = 0.3A and two values of the transmission: (a) 
r = 0.95 and (b) r = 0.8. The different curves correspond to 
different values of a, as indicated in the graphs. The inelastic 
broadening used in these calculations is r\ — 10 _3 A. 



context of atomic contacts^ semiconductor nanowiresi 3 - 
or graphene junctions*^— Therefore, these results in- 
dicate that the microscopic theory presented here will 
always be necessary for the description of the experimen- 
tal results of highly transmissive junctions at sufficiently 
high power, no matter how low the microwave frequency 
is. 



B. Finite temperature: Enhancement of the 
supercurrent 

We now turn to the analysis of the supercurrent at 
finite temperature, carried out within the microscopic 
model. The new ingredient at finite temperature is the 
fact that the ABSs are neither fully occupied nor fully 
empty, which means that quasiparticle transitions be- 
tween the continuum of states and the bound states are 
possible, even for frequencies Huj < A. This has impor- 
tant consequences. 

In Fig. [7] we show the CPR for hu = 0.6A, a = 0.1, 
t = 0.95 and two different temperatures. For compari- 
son, we also show the results in the absence of microwaves 
(dashed lines). Apart from the dips, whose origin is dis- 
cussed above in detail, one can observe that at a certain 
value of the phase ((po ps 0.787r) the current is suppressed. 
Notice that the suppression is stronger as the temper- 
ature increases. Moreover, for phases smaller than <po 
the supercurrent exceeds its value in the absence of mi- 



0.2 0.4 0.6 0.8 1 
Cp/7C 

FIG. 7. (Color online) The current-phase relation for hw = 
0.6A, a = 0.1, t — 0.95 and two different temperatures: (a) 
k B T = 0.2A and (b) k B T = 0.3A. The solid lines in both 
panels correspond to the results of the microscopic theory 
and the dashed lines to the supercurrent in the absence of 
microwaves (a = 0). 



crowaves. In other words, for ip < po there is an enhance- 
ment of the supercurrent induced by the microwave field. 
The origin of this enhancement is the promotion of quasi- 
particles from the continuum below — A to the lower ABS 
by the microwave field. There is also an identical con- 
tribution coming from transitions connecting the upper 
state and the continuum above +A. At low microwave 
powers, these processes can only occur if the field fre- 
quency is larger than the distance in energy between the 
gap edges and the nearest ABS, i.e., if fku > A — Ea(<p), 
and they become possible at finite temperature because 
the lower state is not fully occupied and the upper state 
is not fully empty. For the parameters of Fig. [7] the 
previous condition is satisfied if E\(p) > 0.4A, which 
corresponds to a phase (p < 0.78-7T. Obviously, this phe- 
nomenon of microwave-enhanced supercurrent cannot be 
described by the two-level model since this models ig- 
nores the contribution of the continuum part of the spec- 
trum. 

In order to confirm our interpretation of the origin of 
the microwave-enhanced supercurrent, we have derived 
analytical results describing this phenomenon in the limit 
of weak microwave fields. We have obtained such re- 
sults with the help of an alternative method, known as 
Hamiltonian approach, which for SQPCs has been shown 
to be equivalent to the microscopic theory described at 
the beginning of this section ! 39 ' 42 ' 43 In this approach, a 
point contact is described in terms of a tight-binding-like 
Hamiltonian and the transport properties are calculated 
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following a perturbative approach, where the coupling 
between the electrodes is treated as the perturbation. 
Although the calculations with this method are slightly 
more cumbersome than with the approach used above, 
it has certain advantages. For instance, it also allows 
us to obtain the density of states (DOS) at the contact. 
Moreover, a perturbative analysis (in the field) is much 
simpler when using Eq. (14"5|) . The technical details of the 
Hamiltonian approach are described in the Appendix A, 
and in what follows, we only discuss the results of this 
analysis. 

We are interested in the correction to the current due 
to the microwave field which is responsible for the en- 
hancement of the supercurrent. Thus, based on our nu- 
merical results, we explore the parameter region where 
A — Ea < hjJ. Moreover, in order to avoid the reso- 
nant transitions between the ABSs, we also assume that 
huj < 2Ea- As described in Appendix A, a perturbative 
analysis to lowest order in the parameter a shows that 
the supercurrent can be written as 



-E. 



I(<p)=I«i(<p)+6I(<p), 



(49) 



where I cq is the equilibrium supercurrent given by 
Eq. ([3]), and the correction SI contains several contri- 
butions of order a 2 . There are two types of contribu- 
tions. One type of contribution is related to the change 
in the bound states induced by the coupling to the elec- 
tromagnetic field. The other contribution comes from 
the modification of the occupations of the bound states 
due to the quasiparticle transitions involving the ABSs. 
In the range of parameters that we are interested in, the 
second type of contributions dominates at high enough 
temperatures and, in particular, they are responsible for 
the supercurrent enhancement. Those contributions can 
be written in the spirit of Eq. ([3]) as 



Wenh(v) = Ix(V>) 5n (<?) + / A(^) <5n+ ( ( y 5 )' 



(50) 



where lf(ip) = (2e/h)dE%/dcp give the contribution of 
the states to the equilibrium supercurrent, and Sn^((p) 
are the corrections to the occupations of the ABSs due to 
the application of the microwave field. These corrections 
can be written as 

Sn^ip) = ^ [Re{e i «>p L (E±)p R (E±±oj)} + 

i/(^M^±w)][F -F ±1 ]. (51) 

Here, F n is the distribution function with shifted argu- 
ments F n — ta,nh[(E+nh(jj)/2kBT], v(E) is the density of 
states at the contact in the absence of microwaves, and pj 
and pj are the real part of the anomalous Green functions 
on the left (L) and (R) side of the interface (j = L, R) 
without the field, as defined in Appendix A. Equation 
(|5"Tj) has a very appealing form and it tells us that the 
occupations of the ABSs can be changed by microwave- 
induced transitions connecting these states between the 
continua below and above the gap. These transitions are 
illustrated in Fig. [51 where we also present an example 
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FIG. 8. (Color online) The local density of states at the con- 
tact in the absence of microwaves, as defined in Eq. (|52[) . as 
a function of energy for r = 0.95, <p = 3tt/4 and r\ = 10 _3 A. 
The lower arrows represent the microwave- induced transitions 
between the continuum part of the spectrum and the Andreev 
bound states which are responsible for the supercurrent en- 
hancement at finite temperatures. The upper arrow indicates 
the resonant transition between the ABSs, which suppresses 
the supercurrent. 



of the density of states of the contact in the absence of 
microwaves, v (E). This density of states is given by (see 
Appendix A) 



v(E) = Re 



Ey/{E + ir]) 2 - A 2 

(e + i V y 



E\ 



(52) 



where the poles correspond to the ABSs and, as one can 
see in Fig. [5J there are no singularities at the gap edges 
E = ±A. 

From Eq. (|5"Tj) one can show that the transitions be- 
tween the continuum of states below —A and the lower 
ABS increase the population of the lower state (Sn~ > 0), 
while the photon processes connecting the continuum 
above +A and the upper ABS decrease the occupation 
of the upper state (8n + < 0). As one can see from 
Eq. ([ST))) , both types of processes give a positive contribu- 
tion to the current at finite temperatures and thus, they 
are responsible for the supercurrent enhancement. In- 
deed, due to the electron-hole symmetry of this problem, 
terms in Eq. (|50[) give the same contribution to the cur- 
rent. Finally, the correction to the current due to these 
microwave-induced transitions involving the continuum 
can be written as 



SI e nh((p) = 



-2eE' A 



T 

16 



(53) 



y/(E A + Sw) 2 - A 2 N /A 2 - E 



t]HujEa(2Ea + hu) 
[E A fuj + A 2 (l + cos<^)] [Fi - F ] 6 (\E A + hu\ - A) . 

This expression gives a positive contribution to the su- 
percurrent and it explicitly shows that the enhancement 
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FIG. 9. (Color online) The dc Josephson current as a function 
of the frequency cj of the microwave field for a fixed value of 
the phase if = n/2, and a — 0.1, r = 0.95 and fcsT = 
0.4A. The solid line shows the exact numerical result while 
the dashed line shows the result obtained from Eq. (|53[) . The 
dotted line shows the value of the current in the absence of 
the microwave field. 



can only take place when Hlj > A — Ea- According to 
Eq. (|53p . the correction to the current is proportional 
l/rj, the inelastic scattering time. In our model the pa- 
rameter r\ describes the energy loss mechanism via which 
the microwave power is dissipated. For simplicity, we as- 
sume it to be energy and frequency independent. Equa- 
tion (|53p reproduces the exact results obtained with the 
microscopic approach in the limit of weak fields and in 
the range of frequencies where the transitions between 
the ABSs cannot take place. This is illustrated in Fig. [S] 
where we show the supercurrent for a fixed value of the 
phase (if = 7r/2) as a function of the frequency for 
t = 0.95, a = 0.1 and k^T = 0.4A. As one can see, 
the exact result (solid line) remains constant for low fre- 
quencies. Then, at Hlu = A — Ea there is a rise of the 
supercurrent due to the onset of the transitions connect- 
ing the ABSs with the continuum of states. This increase 
of the supercurrent is well described by the analytical re- 
sult of Eq. (|53|) (dashed line). At higher frequencies, one 
can observe the dips due to the transitions between the 
ABSs. The dip at hw = Ea corresponds to a two-photon 
process, while the one at Huj = 2Ea is produced by a 
single-photon process. Finally, at Huj = A + Ea the su- 
percurrent starts to decrease due to microwave-induced 
transitions between the continuum below — A and the up- 
per ABS and similar ones between the continuum above 
+A and the lower ABS. These transitions, which can also 
occur at zero temperature, tend to increase the occupa- 
tion of the upper state and to reduce the population of 
the lower one, which results in a reduction of the net 
supercurrent. 

As one can see in Fig. [7] and [9j the maximum super- 
current sustained by the junction, i.e. the critical cur- 
rent, can also be enhanced by the microwave field at fi- 
nite temperatures. A microwave-enhanced critical cur- 



rent was first reported in experiments on superconduct- 
ing microbridge o 44 i 45 and explained by Eliashberg 4 ^ in 
1970 in terms of the stimulation of the superconductivity 
in the electrodes, which were made of thin films. Such a 
stimulation, and the corresponding microwave-enhanced 
critical current, only occur at temperatures very close to 
the critical temperature. Enhancements at much lower 
temperatures were reported in the 1970's in the context 
of SNS structures ; 4 ^ 48 and they have been recently ex- 
plained in terms of the redistribution of the quasiparti- 
cles induced by the field42 In this case, for the enhance- 
ment to occur, the temperature must be of the order 
of the minigap in the normal wire, which can be much 
lower than the critical temperature of the superconduct- 
ing leads. 

As discussed above, in the case of a point con- 
tact the mechanism is similar to that of diffusive SNS 
structures^ but it involves discrete ABSs, rather than 
a continuous band of ABSs, as in the case of diffusive 
proximity structures. For this reason we may expect the 
enhancement of the critical current in SQPCs to occur at 
intermediate temperatures, when k^T is of the order of 
the energy distance between the ABSs and the gap edges 
(A — -^A^max)), where ip m ax is the phase value at which 
the supercurrent reaches its maximum. This is illustrated 
in Fig. 1101 where we show the critical current as a func- 
tion of a for different temperatures and different values of 
the transmission. Panel (a) shows the critical current for 
a highly transmissive channel (r = 0.97) and three values 
of the temperature. Notice first that at finite tempera- 
tures, the critical current at finite a (a < 0.5) exceeds the 
value in the absence of microwaves (a = 0). Notice also 
that as a increases, the critical currents clearly deviate 
from the behavior described by the adiabatic approxi- 
mation, which is shown as dashed lines. It is also im- 
portant to emphasize that the microwave-enhancement 
of the critical current is not exclusive of high conductive 
channels and it persists up to relatively low transmis- 
sions, as we show in Fig. HUT b-c). The relative enhance- 
ment of the critical current is larger the larger is the tem- 
perature. It is also worth remarking that at sufficiently 
high power, the critical current depends only weakly on 
the temperature. 



V. CONCLUSIONS AND FURTHER 
DISCUSSIONS 

Summarizing, we have presented a theoretical analy- 
sis of the supercurrent in a phase-biased SQPC under 
microwave irradiation. We have shown that if the mi- 
crowave frequency uj is not high enough to induce tran- 
sitions between the ABSs or between the ABSs and the 
continuum of states outside the gap region, the supercur- 
rent is correctly described by the standard adiabatic ap- 
proximation (see Section II). However, when Huj is com- 
parable to the Andreev gap (energy distance between 
the ABSs), quasiparticle transitions between the ABSs 
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FIG. 10. (Color online) The critical current as a function of 
a for hjj = 0.6 A. The different curves correspond to different 
values of the temperature and the transmission as indicated 
in the graphs. The solid lines correspond to the exact re- 
sults, while the dashed lines show the results of the adiabatic 
approximation. In the three panels the critical current has 
been normalized by eA(T)/h, where A(T) is the gap at the 
corresponding temperature. 



can occur and the supercurrent can be largely suppressed 
at the corresponding values of the phase difference. We 
have shown that this phenomenon can be nicely explained 
within a two-level model that describes the dynamics of 
the ABSsi^i This model indicates that the supercurrent 
suppression is due to the enhancement of the occupation 
of the upper ABS induced by resonant transitions from 
the lower state. Moreover, at low temperatures and weak 
fields, this model is quantitatively correct provided that 
(i) the microwave frequency is not high enough to induce 
transitions connecting the ABSs and the continuum of 
states, and (ii) the Andreev gap is large compared to 
the broadening acquired by the ABSs by means of the 
coupling to the electromagnetic field. Finally, we have 
shown that whenever microwave-induced transitions be- 
tween the ABSs and the continuum of states become pos- 
sible (due to finite temperatures, high frequencies or high 
radiation powers), a fully microscopic theory is needed 



to describe the supercurrent. We have developed such 
a theory and predicted the following effects. First, at 
finite temperatures it is possible to enhance both the su- 
percurrent and the critical current by the application of 
a microwave field. This effect originates from the quasi- 
particle transitions between the ABSs and the contin- 
uum of states, which increase the occupation of the lower 
Andreev state and reduce the population of the upper 
one. Second, the current-phase relation at high pow- 
ers is strongly distorted and it can become highly non- 
sinusoidal exhibiting sign changes in the region between 
and it. Third, the critical current as a function of the ra- 
diation power can exhibit large deviations from the stan- 
dard Bessel-function behavior described by the adiabatic 
approximation. 

It is now pertinent to discuss the connection with ex- 
periments. As explained in the introduction, most of the 
experimental results of the effect of microwaves on the 
supercurrent of a point contact have been successfully 
described in the frame of the adiabatic approximation. 
The reason is that the typical frequency used in the ex- 
periments is relatively low (fioj <C A) and no transitions 
between the ABSs can occur. However, it is important 
to remark that there are no fundamental limitations to 
study the parameter regime where we predict the occur- 
rence of novel effects like the appearance of supercur- 
rent dips in the current-phase relation or the microwave- 
enhanced critical current. These effects are easier to 
observe in highly transmissive point contacts where the 
Andreev gap can become relatively small (much smaller 
than A). The ideal experimental system where to test 
our predictions is a superconducting atomic contact for 
several reasons. First, these contacts can sustain a re- 
duced number of channels, which facilitates the compar- 
ison with the theory. Second, it has been shown that it 
is possible to determine independently the set of trans- 
mission eigenvalues {t^}^ which has allowed to estab- 
lish a comparison between theory and experiment with 
no adjustable parameters for many different transport 



properties 



8,41.51 



. Third, it is possible to tune, at least to 



a certain extent, the transmission coefficients and, in par- 
ticular, to achieve very high transmission coefficients, as 
demonstrated in the context of Al atomic contacts, 8 ' 41 ! 50 
Finally, it has already been shown that in these systems 
the current-phase relation is amenable to measurements, 8 
and investigations of the transport properties of super- 
conducting atomic contacts under microwave irradiation 
have already been performed, 41 ! 52 

In experiments with superconducting atomic junctions, 
even at the level of a single-atom contact, one often 
has the contribution of several conduction channels. In 
this sense, one may wonder whether the presence of low- 
transmissive channels can mask some of the striking ef- 
fects that we have discussed above. In Fig. QT]we show 
the CPR for a contact consisting of three conducting 
channels with transmissions r = 0.17,0.6,0.97 respec- 
tively. The current is obtained by adding the contribu- 
tion of each channel according to Eq. (dJ). As one can 
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FIG. 11. (Color online) The zero-temperature current-phase 
relation for a point-contact consisting of three channels with 
transmissions r = 0.17, 0.6, 0.97, for a = 0.1 and hui = 0.6A. 
The dashed lines show the contribution of each channel, while 
the solid line corresponds to the total current. 
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Appendix A: Hamiltonian approach 

The transport properties of a microwave-irradiated 
SQPC can also be described within the so-called Hamil- 
tonian approachj 42 i 54 We explain in this appendix how 
this approach can be used to obtain analytical results for 
the supercurrent enhancement discussed in Section [IV Bl 
In this approach a single-channel SQPC can be described 
in terms of the following tight-binding-like Hamiltonian 



H = Ht 



H, 



*}, (Al) 



see in Fig. [TT]the total current still shows the dips at the 
resonances corresponding to the channel with the high- 
est transmission (r = 0.97). However, the current does 
not vanish completely due to the contribution of the low- 
transmissive channels. 

It is worth stressing that the major problem to es- 
tablish a direct comparison between our theory and the 
experiments is the fact that we have assumed a phase- 
biased junction. In reality, and depending on the de- 
tails of the electromagnetic environment seen by the 
point contact, the phase across the junction may undergo 
fluctuations (both classical and quantum) which can af- 
fect the value of the critical current or the shape of the 
current-phase relation. Thus, a quantitative comparison 
with the experiments may require in some cases to com- 
bine our theory with a description of the phase fluctua- 
tions. For classical fluctuations, this could be done in the 
spirit of Ref . [53| by means of an extension of the resistively 
shunted junction using our microscopic current-phase re- 
lation as a starting point. 

Let us conclude by saying that in this work we have 
shown that the application of microwaves to one of the 
simplest superconducting systems, namely a SQPC, leads 
to a very rich phenomenology, which has remained largely 
unexplored. In particular, we have shown that a mi- 
crowave field is an ideal tool to make a direct spec- 
troscopy of the Andreev bound states of a superconduct- 
ing junction. The ideas put forward in this work pave the 
way for the understanding of the influence of a microwave 
field on the supercurrent of variety of highly transmissive 
superconducting weak links. 



where Hl,r are the BCS Hamiltonians describing the left 
(L) and right (R) electrodes and the last term describes 
the coupling between the electrodes. In this last term, t 
is a hopping element that determines the transmission of 
the contact. 

In this model the current evaluated at the interface 
between the two electrodes adopts the form 



^) = |E{*< c L^ ff }-r<cL^>} 



(A2) 



This expression can be rewritten in terms of the Keldysh 
Green functions as 



(A3) 



Here T3 is the corresponding Pauli matrix, Tr denotes 
the trace in Nambu space, and t is the hopping matrix 
in Nambu space given by 



t = 



ie#(t)/2 






_ t * e -i4>(t)/2 



(A4) 



Here, <f>(t) is the time-dependent superconducting phase 
given by Eq. ([2]) . 

In order to determine the Green functions appearing in 
the current expression, we follow a perturbative scheme 
and treat the coupling term in Hamiltonian (|A1|) as a 
perturbation. The unperturbed Green functions describe 
the uncoupled electrodes in equilibrium. Thus for in- 
stance, the retarded and advanced functions are given 
by 



g n {*>) - 



i 



WC R(A) (E) 



E A 
A E 



(A5) 



15 



where j = L,R, C R(A) = ^(E + irj) 2 - A 2 , and W is 
an energy scale related to the normal density of states 
at the Fermi energy. The full Green functions can then 
be determined by solving a Dyson equation, where the 
retarded and advanced self-energies are simply given by 
the hopping matrix of Eq. (|A4[) . 

Since we are interested in the limit of weak fields (a <C 
1), we can expand the phase factors in Eq. (|A4[) as follows 



J<t>(t)/z 



1 + a cos ujt 



1 



-a 2 (cosu;i) 2 



Moreover, for the perturbative treatment in a it is conve- 
nient to use the full Green functions of the contact in the 
absence of microwaves (a = 0), Gij. It is straightforward 
to show that these functions can be expressed as 



G 



R{A) 
LL 



W(l + (3)Z R ( A ) V E a E±i v 



E±ir) 



e; 



(A7) 



A,R{A) = -t 

W 2 (l + j3)£ R W 



a R(A) 
_ h R{A)* 



h R(A) 
_qR(A)* ) 



where E g = A(l+^)(l+/3), /3 = (t/W) 2 , £ = E 2 -E\, 
a = E 2 e- i '*'/ 2 -AE*e i ^ 2 and b = E{E g e-^/ 2 - Ae iip ^). 
Similar expressions hold for Grr and Glr- These Green 
functions are now the zero-order propagators of the per- 
turbation theory. Substituting these functions in the cur- 
rent expression of Eq. (|A3[) and identifying the trans- 
mission coefficient as r = 4/3/(1 + /3) 2 ,— one obtains 



the expression for the equilibrium current of Eq. ([3]). 
On the other hand, from the previous expressions one 
can determine the local density of states at the con- 
tact in the absence of microwaves, which is defined as 
Vj (E) = W(l + /9)(i/2)(6£. - 6%)u {vl = vr in our 
symmetric contacts). This density of states is given by 
Eq. (|52|) and it is shown in Fig.! 

Going into the energy representation as done in Section 
TVl the first correction to the current of Eq. (|A3[) . which 
is of order a 2 , contains the following three terms 



51 = 



dE 
2^ 



where the superindices denote the order of perturbation 
in a. To obtain a complete analytical expression for an 
arbitrary value of the field frequency is quite cumber- 
some. Instead, we concentrate on the parameter range 
where the current enhancement takes place. For that 
purpose, we focus on frequency values far from the reso- 
nant condition tuo = 2Ea and close to A — Ea- In this 
region, it turns out that the second term in Eq. (|A9|) is 
proportional to the parameter A/?/. All the other terms 
give a contribution which depends only weakly on the 
frequency. Assuming a small inelastic scattering rate, 
one can approximate the correction to the current by 
Eqs. & and HE}, with Pj = (i/2)W(l + f3)[Gf-Gf] h2 , 

Pj = {i/2)W(l + p)[G?-Gf} 2 ,i , andF n = F(E + nhw). 
This correction gives precisely the enhancement of the su- 
percurrent at finite temperatures, as discussed in section 
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